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j> ! Abstract: We consider low scale slow roll inflation driven by the gauge invariant 

•1— I , 

^ ■ flat directions udd and LLe of the Minimally Supersymmetric Standard Model at 

the vicinity of a saddle point of the scalar potential. We study the stability of saddle 



point and the slow roll regime by considering radiative and supergravity corrections. 
The latter are found to be harmless, but the former require a modest finetuning of 
the saddle point condition. We show that while the inflaton decays almost instantly, 
full thermalization occurs late, typically at a temperature T w 10 7 GeV, so that 
there is no gravitino problem. We also compute the renormalization group running 
of the inflaton mass and relate it to slepton masses that may be within the reach of 
LHC and could be precisely determined in a future Linear Collider experiment. 
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1. Introduction 



Recently we have constructed a model of inflation |l[ based on the udd and LLe 
flat directions of Minimally Supersymmetric Standard Model (MSSM; for a review 
of MSSM flat directions, see f2|). In this model the inflaton is a gauge invariant 
combination of either squark or slepton fields. For a choice of the soft SUSY breaking 
parameters A and the inflaton mass m^, the potential along the flat udd and LLe 
directions is such that there is a period of slow roll inflation of sufficient duration to 
provide the observed spectrum of CMB perturbations. In the inflationary part of the 
MSSM potential the second derivative is vanishing and the slow roll phase is driven 
by the third derivative of the potential 1 . 

MSSM inflation occurs at a very low scale with H in f ~ 1 — 10 GeV and with 
field values much below the Planck scale. Hence it stands in strong contrast to the 
conventional inflation models which are based on ad hoc gauge singlet fields and 
often employ field values close to Planck scale (for a review, see Q). In such models 
the inflaton couplings to SM physics are unknown. As a consequence, much of the 
post-inflationary evolution, such as reheating, thermalization, generation of baryon 
asymmetry and cold dark matter, which all depend very crucially on how the inflaton 
couples to the (MS)SM sector [|, || |7|, §, is not calculable from first principles. The 
great virtue of MSSM inflation based on flat directions is that the inflaton couplings 
to Standard Model particles are known and, at least in principle, measurable in 
laboratory experiments such as LHC or a future Linear Collider. 

However, as in almost all inflationary models, a fine tuning of the initial condition 
is needed to place the flat direction field to the immediate vicinity of the saddle 
point <fto at the onset of inflation. In addition, there is the question of the stability 
of the saddle point solution and of the existence of a slow roll regime. These are 
issues that we wish to address in detail in the present paper. Both supergravity 
and radiative corrections to the flat direction inflaton potential must be considered. 
Hence we need to write down and solve the renormalization group (RG) equations 
for the MSSM flat directions of interest. RG equations are also needed to scale the 
model parameters, such as the inflaton mass, down to TeV scale; since the inflaton 
mass is related either to squark or slepton masses, it could be measured by LHC or 

a recent similar model with small Dirac neutrino masses, the observed microwave background 
anisotropy and the tilted power spectrum are related to the neutrino masses ||. The model relies 
solely on renormalizable couplings. 



-2- 



a future Linear Collider. 

Because the inflaton couplings to ordinary matter are known, inflaton decay and 
thermalization are processes that can be computed in an unambiguous way. Unlike 
in many models with a singlet inflaton, in MSSM inflation the potential relevant for 
decay and thermalization cannot be adjusted independently of the slow roll part of 
the potential. 

This paper is organized as follows. In Sect. 2 we present the model of MSSM 
inflation and its predictions. In Sect. 3 we study the flat direction potential without 
an exact saddle point. We find generic constraints for the existence of a slow roll 
solution and show that in the slow roll regime there is always tunneling from a false 
minimum. In Sect. 4 we solve the one-loop RG equations for the LLe flat direction 
and find the one- loop corrected saddle point. We quantify the amount of fine tuning 
required for the slow roll solution to exist, and relate through RG running the LLe 
inflaton mass with observables such as the slepton masses at the LHC energy scale. 
We also show that supergravity corrections to the potential can be neglected. In 
Sect. 5 we discuss the decay of the flat direction, the reheating and thermalization 
of the Universe, and show that the reheat temperature is low enough for the model 
to avoid the gravitino problem. Sect. 6 contains our conclusions and some discussion 
about future prospects. 



2. The Model 

Let us recapitulate the main features of MSSM flat direction inflation [I]]. As is 
well known, in the limit of unbroken SUSY the flat directions have exactly vanishing 
potential. This situation changes if we take into account soft SUSY breaking and 
non-renormalizable superpotential terms 2 of the type @ 

n>3 

where $ is a superfield which contains the flat direction. Within MSSM all the flat 
directions are lifted by non-renormalizable operators with 4 < n < 9 [|TTJ, where n 



2 Our framework is MSSM together with gravity, so consistency dictates that all non- 
renormalizable terms allowed by gauge symmetry and supersymmetry should be included below 
the cut-off scale, which we take to be the Planck scale. Some interseting issues on A-term inflation 
were also discussed in Ref. IrUl. 
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Figure 1: The colored curves depict the full potential, where V(x) = 
Vf»/(0.5 mlM^m^/Mp) 1 / 2 ), and x = (A n M P /m (} i) 1 / 4 (0/Mp). The black curve is the 
potential arising from the soft SUSY breaking mass term. The black dots on the colored 
potentials illustrate the gradual transition from minimum to the saddle point and to the 
maximum. 

depends on the flat direction. We expect that quantum gravity effects yield M = 
M P = 2.4 x 10 18 GeV and A n ~ 0(1) 3 . 

Let us focus on the lowest order superpotential term in Eq. (|2.1|) which lifts the 
flat direction. Soft SUSY breaking induces a mass term for and an A-term so that 
the scalar potential along the flat direction reads 

V = -ml^ + A cos(ne + 9 A ) + A* , (2.2) 

Here <ft and 9 denote respectively the radial and the angular coordinates of the com- 
plex scalar field $ = exp[z#], while 9 a is the phase of the A-term (thus A is a 
positive quantity with dimension of mass). Note that the first and third terms in 
Eq. Q2.2p are positive definite, while the A-term leads to a negative contribution 
along the directions whenever cos(n9 + 9a) < 0. 4 

3 Note however that our results will be valid for any values of A„, because rescaling A„ simply 
shifts the VEV of the flat direction. 

4 The importance of the A-term has also been highlighted in a successful MSSM curvaton 
model || . 
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In principle, in the A-term all the superpotential terms of a given dimension n 
may enter with a different coefficient A n ; whether they are related or not depends 
on the details of the SUSY breaking mechanism. 

2.1 The Saddle Point 

The maximum impact from the A-term is obtained when cos(n6 + 9a) = —1 (which 
occurs for n values of 9). Along these directions V has a secondary minimum at 
= O ~ (m^Mp -3 ) 1 ^ 2 <C Mp (the global minimum is at = 0), provided that 

A 2 > 8(n - l)m^ . (2.3) 

At this minimum the curvature of the potential is positive both along the radial and 
angular directions 5 with V ~ ^ 2 0o ~ w| (m^Mp" 3 ) 2 "" 2 \ 

As discussed in [|IJ, if the local minimum is too steep, the field will become 
trapped there with an ensuing inflation that has no graceful exit like in the old 
inflation scenario |TJ]]. On the other hand in an opposite limit, with a point of 



inflection, a single flat direction cannot support inflation | 15fl . 

However, in the gravity mediated SUSY breaking case, the A-term and the soft 
SUSY breaking mass terms are expected to be the same order of magnitude as the 
gravitino mass, i.e. 

~ A ~ m 3/2 ~ 0(1) TeV . (2.4) 

Therefore, as pointed out in |TJ, in the gravity mediated SUSY breaking it is possi- 
ble that the potential barrier actually disappears and the inequality in Eq. (|2.3|) is 
saturated so that A and are related by 

A 2 = 8(n-l)mJ. (2.5) 

This represents a fine tuning and will be discussed at length in the next Sections. 
However, let us now assume for the sake of argument that Eq. ( |2.5| ) holds. Then 
both the first and second derivatives of V vanish at O , i.e. ^'(^o) = 0, V"(0o) = CI- 
As the result, if initially ~ 0o, a slow roll phase of inflation is driven by the third 
derivative of the potential. 



5 If the A is too large, the secondary minimum will be deeper than the one in the origin, and 
hence becomes the true minimum. However, this is phenomenologically unacceptable as such a 
minimum will break charge and/or color fll2fl. 
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Note that this behavior does not seem possible for other SUSY breaking scenarios 
such as the gauge mediated breaking |L6) or split SUSY [fT7|| . In split SUSY the A- 



term is protected by an i?-symmetry, which also keeps the gauginos light while the 
sfermions are quite heavy |T7| 6 . 



2.2 Slow roll 

The potential near the saddle point Eq. ( |2.5| ) is very flat along the real direction but 
not along the imaginary direction. Along the imaginary direction the curvature is 
determined by m^. Around 0o the field lies in a plateau with a potential energy 

with 



X n y/2n - 2, 

This results in Hubble expansion rate during inflation which is given by 



(2.7) 



(n - 2) m^o 
inf ~ y/6n(n - 1) M P " ^ 

When is very close to <po, the first derivative is extremely small. The field is 
effectively in a de Sitter background, and we are in self-reproduction (or eternal 
inflation) regime where the two point correlation function for the flat direction fluc- 
tuation grows with time. But eventually classical friction wins and slow roll begins 
at (p « se i f 

(0o-0sci f )-(^) 1/2 0o. (2.9) 
The slow roll potential in this case reads 

V{<f>) = V(<h) + 7^'"(0o)(0 - 0o) 3 + • ■ ■ , 

9 

V /,/ (0o) = 2(n-2) 2 -/. (2.10) 
We can now solve the equation of motion for the field in the slow-roll approximation, 

3i?0=-iy'"(0 o )(0-0o) 2 , (2.11) 



3 In the gauge mediated case there is an inherent mismatch between A and m^, except at very 



large field values where Eq. (2.4) can be satisfied. However there exists an unique possibility of a 



saddle point inflation which we will discuss separately J18[ 
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assuming initial conditions such that the flat direction starts in the vicinity of 0o with 
rs 0. Inflation ends when the slow roll parameter, e = (Mp/2)(V /V) 2 becomes of 
(9(1). This occurs at 



;,3 

~ — TTT72 • (2-12) 



2n(n - l)Mj 

which happens to be also the place when the other slow roll paremeter rj = Mp(V" /V) 
becomes of 0(1). 

The number of e-foldings during the slow roll from <p to en d is given by 

KW =r^ ., JLn it ' (2.13) 



2n(n - 1)M|(0 O - 

where we have used V(</>) ~ (0 — (fio) 2 V" (4>o) (this is justified since V'{4>q) ~ 
0, V"(0o) ~ 0), and Eq. ( |2.11| ). The total number of e-foldings in the slow roll 
regime is then found from Eq. (|2.9|) 

Mot-^— pr-^r) • (2-14) 

The observationally relevant perturbations are generated when rs 0cobe- The num- 
ber of e-foldings between 0cobe and 4> en d, denoted by Acobe follows from Eq. ( ^.1^ ) 

^COBE ^ 7T7 TTtJt 7 V- ( 2 - 15 ) 

2n(n - l)M^(0o - </>cobe) 
The amplitude of perturbations thus produced is given by 



where we have used Eqs.( [2.8|) , (|2.10|) , ( |2.15|) . Again after using these equations, the 
spectral tilt of the power spectrum and its running are found to be 

4 

n s = 1 + 2r] - 6e ~ 1 - — : , (2.17) 

•A/cOBE 

dn s 4 , 

(2.18) 



dink A^ OBE ' 
2.3 Properties and predictions 

As discussed in [[|], among the about 300 flat directions there are two that can lead 
to a successful inflation along the lines discussed above. 

One is udd which, up to an overall phase factor, is parameterized by 

< = i=0, d?=-^, 4=^- (2-19) 
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Here 1 < at, /3, 7 < 3 are color indices, and 1 < i, j, k < 3 denote the quark families. 
The flatness constraints require that a 7^ (3 7^ 7 and j 7^ A;. 

The other direction is LLe, parameterized by (again up to an overall phase 
factor) 

> (2 - 2o> 

where 1 < a, 6 < 2 are the weak isospin indices and 1 < i, j, A; < 3 denote the lepton 
families. The flatness constraints require that a 7^ b and i ^ j ^ k. Both these flat 
directions are lifted by n = 6 non-renormalizable operators, 

W 6 D -^(LLe)(LLe) , W 6 D -^(udd)(udd) . (2.21) 

The reason for choosing either of these two flat directions 7 is twofold: (i) a non- 
trivial A-term arises, at the lowest order, only at n = 6; and (ii) we wish to obtain 
the correct COBE normalization of the CMB spectrum. 

Those MSSM flat directions which are lifted by operators with dimension n = 7, 9 
are such that the superpotential term contains at least two monomials, i.e. is of the 
type 

W ^^$ n - 1 . (2.22) 

Mp V ' 

If represents the flat direction, then its VEV induces a large effective mass term for 
if), through Yukawa couplings, so that (■0) = 0. Hence Eq. ( |2.22| ) does not contribute 
to the A-term. 

More importantly, as we will see, all other flat directions except those lifted by 
n = 6 fail to yield the right amplitude for the density perturbations. Indeed, as can 
be seen in Eq. ( |2.7] ), the value of 0o, and hence also the energy density, depend on n. 

According to the arguments presented above, successful MSSM flat direction 
inflation has the following model parameters: 



lOTeV, n = 6, A = V^Om^ , A ~ 0(1) . (2.23) 



Here we assume that A (we drop the subscript "6") is of order one, which is the most 
natural assumption when M = Mp. 



7 Since LLe are udd are independently D- and F-flat, inflation could take place along any of 
them but also, at least in principle, simultaneously. The dynamics of multiple flat directions are 
however quite involved JL9[. 
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The Hubble expansion rate during inflation and the VEV of the saddle point 
are 8 

H in{ -1-10 GeV , 0o ~ (1 - 3) x 10 14 GeV . (2.24) 

Note that both the scales are sub-Planckian. The total energy density stored in the 
inflaton potential is Vq ~ 10 36 — 10 38 GeV 4 . The fact that 4>q is sub-Planckian guar- 
antees that the inflationary potential is free from the uncertainties about physics at 
super-Planckian VEVs. The total number of e-foldings during the slow roll evolution 
is large enough to dilute any dangerous relic away, see Eq. ( [2.14| ): 

Mot ~ 10 3 , (2.25) 



Domains which are initially closer than se if to (fio, see Eq. ( |2.9| ), can enter self- 
reproduction in eternal inflation, with no observable consequences. 

At such low scales as in MSSM inflation the number of e-foldings, Acobe, re- 



quired for the observationally relevant perturbations, is much less than 60 [2D]. If 
the inflaton decays immediately after the end of inflation, we obtain Acobe ~ 50. 
Despite the low scale, the flat direction can generate adequate density perturbations 
as required to explain the COBE normalization. This is due to the extreme flatness 
of the potential (recall that V = 0), which causes the velocity of the rolling flat 



direction to be extremely small. From Eq. (|2.16|) we find an amplitude of 



5 H ^ 1.91 x 10~ & . (2.26) 

There is a constraint on the mass of the flat direction from the amplitude of the 
CMB anisotropy: 

m ~ (100 GeV) x A" 1 (^p) ' • (2.27) 

We get a lower limit on the mass parameter when A < 1. For smaller values of 
A< 1, the mass of the flat direction must be larger. Note that the above bound on 
the inflaton mass arises at high scales, i.e. = 4>q. However, through renormalization 
group flow, it is connected to the low scale mass, as will be discussed in Sect. 4. 

The spectral tilt of the power spectrum is not negligible because, although the 
first slow roll parameter is e ~ 1 / Acobe ^ 1 > the other slow roll parameter is given 



8 We note that f/i n f and 4>o depend very mildly on A as they are both oc A 
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by f] = — 2/A/cobe and thus, see Eq. (|2.17 ) 



9 



n s ~0.92, (2.28) 

—4 0.002, 2.29 

dmk 

where we have taken A/cobe ~ 50 (which is the maximum value allowed for the scale 
of inflation in our model). In the absence of tensor modes, this agrees with the current 



WMAP 3-years' data within 2a ||22|| . Note that MSSM inflation does not produce 
any large stochastic gravitational wave background during inflation. Gravity waves 
depend on the Hubble expansion rate, and in our case the energy density stored in 
MSSM inflation is very small. 

3. Sensitivity of the saddle point inflation 

In previous Sections and in Ref. [IJ the dynamics of the flat direction inflaton was 
discussed assuming the saddle point condition Eq. (|2~5|) is satisfied exactly. The 
question then is, how large a deviation can be allowed for before slow roll inflation will 



be spoiled. There are obviously two distinct possibilities: either A > \/9>{n — 1)??^ 
or A < \fS(n — l)m0. (Although we always take n = 6 in the present paper, we keep 
n here for generality of the formalism.) In the former case there is a barrier which 
separates the global minimum = and the false minimum at <fi ~ <J)q. The eventual 
inflationary trajectory starts near the top of the barrier. The field can either start 
at the top, or jump to its vicinity from the false minimum via Coleman-de Luccia 
tunneling |23]. As we will see, if the barrier is too high, there will be no inflation 



near the top. In the latter case there is no minimum but the potential may be too 
steep for slow roll inflation. Therefore we need to analyze the two cases separately. 
However, the steepness of the potential is a problem which is common to both cases 
and is addressed at the end of this Section. 
To facilitate the discussion, let us define 

6 = , A \ „ = 1 ± ( —) 2 a 2 . (3.1) 
8(n-l)mJ V 2 J y } 

Here we will assume that a < 1. Before beginning the calculations, we would like to 
point that the main results of this section are summarized in Fig. (|3]) and Eq. ( |3.31| ). 

9 Obtaining n s > 0.92 (or n s < 0.92, which is however outside the 2a allowed region) requires 



deviation from the saddle point condition in Eq. (2.5), see Section 3. For a more detailed discussion 
on the spectral tilt, see also Refs. fl,pi|. 
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These yield no constraint on the spectral tilt as any value consistent with sufficient 
slow roll inflation (i.e. a number > A/cobe of e-foldings) is allowed. 

3.1 The potential for 5 > 1 

In this case there are two extrema, a maximum (— ) and a minimum (+), 



>± = do \V5 ± V5^T] n ' 2 = 0o (l ± ~ + 0(a 2 )) , (3.2) 



l + ^-^a 2 + 0(a 4 )). (3.3) 



and a point of inflection 



32 

We can then express the potential and its derivatives at the extrema 0_ and d> + as 
functions of a, 

V(<f> ± ) = Vfa) (l - ± ^_}l a ^ j (3 . 4) 

K'(0 ± ) = ±a(n-2) 2 mJ, (3.5) 
K"(0 ± ) = K"(0 O ) (l ± ^ - 2)a) . (3.6) 

Note that a — > when 5 — > 1, so we can expand the potential around the maximum, 
and include the small a correction due to deviations from the saddle point as 



V((/>) = (n- 2) 2 m 2 2 



1 a f d> a\ 2 1 / 6 a x ' 



2n(ra- 1) 2 V0o 2/ 3 V0o 2 



(3.7) 

The maximum is now at max = </>o(l — a/2), and the minimum is at </> min = d>o(l + 
a/2), with masses (curvature of the potential) given by m 2 ± = ±(n — 2) 2 m| a, which 
coincides with Eq. ( |3.5| ), while there is now a point of inflection at d>o- Note that the 
difference in potential height between the maximum and the minimum is 

Kiax - Knin 7l(n - 1) 3 
= ; : = Oi . (O.o) 

In the limit of a — > we recover the saddle point. We will work in the limit when 
a < 1. 

Let us now define a few variables, x = d>/d>o, T — (n — 2)m^t and h = 
H- m f/(n — 2)m0. Then the equation of motion for the scalar field down the potential 
can be written as 

x" + 3hx' + V\x) = x" + 3hx' + (x - 1 + |) 2 - a(x - 1 + |) = , (3.9) 
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where we have used 

The eventual inflationary trajectory will start in the vicinity of the maximum 
x = 1 — a/2 and will roll down the hill towards x = 0. The field can either start near 
the maximum, or tunnel to its vicinity out of the false vacuum. Tunneling takes place 
in the presence of a non-zero vacuum energy, H m f ^ 0, and is known as Coleman-de 



Luccia tunneling [23 



In order to find the interpolating solution between the false and the true vacuum 
one solves the Euclidean equation of motion, 

x" + 3hx' -V'(x) = x" + 3hx' - (x - 1 + |) 2 + a(x - 1 + |) = , (3.10) 

whose exact solution is 

Ct CtT 

x(t) = 1 tanh — . (3.11) 

v ' 2 6h K ' 

This solution starts at x(— oo) = 1 + a/2 = (p + /4>o and ends at x(oo) = 1 — a/2 = 
0_/0o- The "tunneling" from min to max can actually be understood as diffusion 
due to de Sitter fluctuations. It is valid so long as V"(0min) "C Hf ni . This requires 
that, see Eq. (|3.5|), 



1 



l2 



a < - / (3.12) 
Qn(n - 1) M| v ; 

This also insures that \t)\ 1 at the maximum, and hence inflation can take place 

after tunneling. Otherwise there will be no inflation, neither in self-reproduction nor 

in slow roll regime. 

Let us now discuss the effect of the tunneling solutions on the tilt of the CMB 
spectrum. Again there is self-reproduction close to the maximum as long as the cur- 
vature there is smaller than the rate of expansion squared, i.e. a -C l/18O(0o/^p) 2 - 
The slow-roll regime starts at 4> ~ 4>seif when 

|(0 self -0_)(0 self -0 + )|<!^o, (3.13) 



Note that |(0 self - <p_)\ ~ |(</> + - <j>_)\ = a<p , see Eq. (|OD, and |(0 se i f - 
|(0seif — 4>o)\- We therefore find 



(3.14) 



Now we integrate Eq. (|3.9| ) in the slow-roll approximation, using a new variable 
u = x — x_, for which the equation of motion becomes 3hu' = —u(u — a). The exact 

solution is, in terms of the number of e-folds, M = hr, 

a 

X (N) = 1 - - cothOaAf) (3.15) 
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Figure 2: The red curve depicts the potential, V(x) = V^/^.hrn^M^m^/M^) 1 / 2 ), 
where x = (X n M-p/m ( j ) ) 1 ^(4>/Mp), for the saddle point (shown by the black dot) when 
5 = 1. The blue curve illustrates the potential when d = l + \/40/1000, where the two black 
dots, one on right shows the minimum value 0_ and on the left shows the maximum value, 
The green curve portrays the potential for the opposite case when 5 = 1 — \/40/1000. 
The black dot is the point of inflection. 



where we have defined 

fi .JL.« n - 1)a !g. (3.16) 

Note that in the limit a — > 0, we recover the usual expression Eq. ( |2.13| ). From 
Eqs. ( |2.12| ), ( |3.14| ), ( |3.15| ), it turns out that the number of e- folds from self to the 
end of inflation at en d is again of order 10 3 . 

The required number of e-folds for the relevant perturbations (Acobe ~ 50) 
determines the value of 0cobe, 

\(M A \(A A Ml/2 a 3/i Vo 00 

0cobe - 0- 0cobe - 0+ 7 = . r o K r - 77 = ^T7 i\m2 K r • 

2 smh fJN Acobe 2n(n-l)M§ Acobe 

(3.17) 

On the other hand, the amplitude of fluctuations is given by 

3 (n-2)^i(2n-2)— ' ( m^V^ sinh 2 /W 

^-s7F(0) = W5 n UrJ — a*" { ] 
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which for n = 6,4,3 becomes 

&H{n = Q) ~ 2 

5 H (n = 4) ~ A n 



A n m^\ 1//2 sinl^/W 

sinh 2 (3X 

J 2 ' 



2 M P sinh 2 /W 



(3.19) 



fc(n = 3)*0.5A^ p2 
while the spectral tilt and its running are universal, 



dln62n l-4/?coth/W, 



dM 

dn q 



4{3 2 



(3.20) 
(3.21) 



dink dM sinh 2 /W 

which reduce to the usual expressions in the limit f3 — > 0, see Eqs. (|2.17|) , ( |2.18| ). 
We show in Fig. |] the variation of the tilt with (3 for a model with JVcobe — 50. 
Note that the range of allowed values of (3 is constrained by the condition to have 
inflation near the maximum, i.e. that \q\ < 1 at 0_. This gives (3 < 0.06, for n = 6, 
see Eqs. ( P . 1 2| , |3 . 1 6|) . The corresponding range of tilt values agrees with the results 
of Ref. §. 

3.2 The potential for 8 < 1 

When 5 < 1, instead of a saddle point we have a point of inflection at <fii, where 
V"(d>i) = 0. We find 



n-2 
32 



a 2 + 0{a 4 



and 



V{<j>i) = V(cf>o) ( 1 + 

2 



3(n-2)(n-4) 
16n 



a 2 + 0{a A 



V'i 



n-2 



a 2 m^ + 0(a 4 



V'"^) = V'Xfa) (1 + 0(a 2 )) . 



The slow-roll parameters at the point of inflection are 



1M 2 fn(n - r 
2 0f 



= 



(3.22) 

(3.23) 

(3.24) 
(3.25) 

(3.26) 
(3.27) 
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Figure 3: The tilt of the MSSM model as a function of the saddle point deviation param- 
eter j3 (3.16). We have plotted both cases 5 < 1 and 5 > 1. Note the allowed 2a range 
coming from a combination of CMB and LSS observations. 

Unlike the previous cases, 5 = 1 (saddle point), and 5 > 1 (tunneling solution), 
there is no point (except the origin = 0) where V = 0. This implies that there 
will be no self-reproduction regime unless a <C (m^/Mp) 1 / 2 |J. 

However this is not troublesome as long as we have a sufficient number of e- 
foldings, Af > Acobe arising due to a slow roll inflation. 

The amplitude and tilt of the scalar spectrum in the case 5 < 1 can be obtained 
from the analytical continuation (3 — > i(3 of the results of previous subsection (5 > 1), 

^3 H*: (n-2)y^(2n-2)^ ^_ ( m<j> \^ s j n 2 (3Af 



n s = 1-4(3 cot ^JV , (3.29) 

-^=- = 4/32 (3 30) 

din A; sin 2 /W' 1 j 

which is in agreement with the results of Ref. 0. 

The dependence of the tilt on j3 can be seen in Fig. 0. Note that, as pointed out 
in Ref. B, the tilt can get any value in the allowed range of (3 which is determined 
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by the viability of slow roll. In the future we will have to determine what value of (3 
agrees with observations. 

To summarize the fine tuning issue, for typical values of 0o ^ 10 15 GeV 10 , the 
saddle point condition, Eq. (|2.5|), requires fine-tuning at the level of 

a~l(T 9 -i(r 8 , (3.31) 

which is not negligible. 



4. Radiative and supergravity corrections 

The MSSM inflaton candidates are represented by gauge invariant combinations but 
are not singlets. The inflaton parameters receive corrections from gauge interac- 
tions which, unlike in models with a gauge singlet inflaton, can be computed in a 
straightforward way. Quantum corrections result in a logarithmic running of the soft 
supersymmetry breaking parameters and A. One might then worry about their 
impact on Eq. ( |2.5|) and the success of inflation. 

In this section we will discuss running of the potential with VEV-dependent 
values of m^cp) and A((p) in Eq. ( |2.5| ). Our conclusion is that the running of the 
gauge couplings do not spoil the existence of a saddle point. However the VEV 
of the saddle point is now displaced; by how much will depend precisely on the 
inflaton candidate. In order to discuss the situation, we derive a general expression 
for the one-loop effective potential for the flat directions, and then focus on the LLe 
direction, for which the system of RG equations can be solved analytically. 

4.1 One-loop effective potential 

The first thing to check is whether the radiative corrections remove the saddle point 
altogether. The object of interest is the effective potential at the phase minimum 



n ^mm — 7T, for which we obtain 



1 

V e ff(<j>,0 min ) = -mltf 



1 + K x log , , 



\2 

, A nfi i2(n-l) 



1 + ^slog 



nM n ~ 3 

2\ 1 



/'0 



1+K 2 log, , 



(4.1) 



where m , A Q , and A nj o are the values of m^, A and A n given at a scale //o- Here A Q 
is chosen to be real and positive (this can always be done by re-parameterizing the 



10 



The tendency from radiative corrections is to raise <^>oj see Section 4.3. 
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phase of the complex scalar field <ft), and \Ki\ < 1 are coefficients determined by the 
one-loop renormalization group equations. 

Our aim is to find a saddle point of this effective potential, so we calculate the 1st 
and 2nd derivatives of the potential and set them to zero. This is a straightforward 
although somewhat cumbersome exercise that results in the expression 



^ 2 



M 



n—3 



iXJn -I + K3 



A ( 1 + ~K 2 J ± J A 2 (l + ^K 2 ) cS//r ( I /v 1 )( // I Is, , 



(4.2) 

where m^, A, and X n are values of the parameters at the scale <pQ. Inserting this into 
y,<j>4> — 0, we can then find the condition to have a saddle point at (f> : 



A 2 = 2m 2 Jn-l + K 3 )F 1 F 2 F 3 



Fx 
F 2 
F 3 



1 + Ki 

n-l + K. 



■(jji - l)(2n - 3) + {in - h)K 3 \ -\-ZKi 



2n — 1 \ (2 

1 + K x ) [ n - 1 + 2 K 2 - (1 + 3^) [I + -K2 

n I \ n 



l + ^K 2 

n L 

n-l + K, 



■(jn - l)(2n - 3) + (in - 5)if 3 ) - f 



2n- 1 
n - 1 + 2 K 2 

n 



-, -1 



(4.3) 



In the limit when \K{\ <C 1, this mercifully simplifies to 



A 2 = 8(n - l)mj(^o) f 1 + i^i - -^2 + -^-r^ 
A rtV ^(^T)V 2 1 2(n-l) 



(4.4) 
(4.5) 



Note that Eqs. (|4.3|) , (|4.4j) give the necessary relations between the values of and 
^4 as calculated at the saddle point. The coefficients Ki need to be solved from the 
renormalization group equations at the scale given by the saddle point u = 0o- Since 
Ki are already one loop corrections, taking the tree-level saddle point value as the 
renormalization scale is sufficient. 

Hence we may conclude that, although the soft terms and the value of the saddle 
point are all affected by radiative corrections, they do not remove the saddle point 
nor shift it to unreasonable values. The existence of a saddle point is thus insensitive 
to radiative corrections. 
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4.2 RG equations for the LLe direction 

The form of the relevant RG equations depend on the flat direction. RG equations 
for LLe are simpler since only the SU(2)w x U(l)y gauge interactions are involved 
and the lepton Yukawa couplings are negligible. The case of udd requires numerics 
if u is chosen from the third family. For other choices, however, it closely resembles 
LLe. For LLe the one-loop RG equations governing the running of rr&, A, and A 
with the scale \i are given by 



dm\ _ 


1 


dfi 


Qtt 2 


dA 


1 


^ dfi 


2tt 2 


dX 


1 


^ dfi 


4tt 2 



3 M o o 3 „, o > 

v—r^ = ~— 1 2 m ^ 2 + 2 mi91 



2 m *92 + ^ m i9i 
*(§«5 + frf)- («) 

Here rhi, m 2 denote the mass of the U(l)y and SU{2) W gauginos respectively and 
gi, #2 are the associated gauge couplings. It is a straightforward exercise to obtain 
the equations that govern the running of A and A associated with the (LLe) 2 su- 
perpotential term (which lifts the LLe flat direction). Note that L has the same 
quantum numbers as ILj, and hence in this respect LLe combination behaves just 
like H^Le. One can then use the familiar RG equations that govern the Yukawa 
coupling and yl-term associated with the H^Le superpotential term |24}| - However, 



as explained in ||25|| , the coefficients of the terms on the right-hand side are propor- 
tional to the number of superfields contained in a superpotential term. 11 Hence the 
second and third equations in ( |4.6| ) are simply obtained from those for the H^Le term 
after multiplying by a factor of 2. The first equation in ( |4.(j| ) is also easily found by 
taking the electroweak charges of Li, Lj and e superfields into account while taking 
into account that m 2 = (yr? L . + m\. + m^)/3. 

The running of gauge couplings and gaugino masses obey the usual equations [EJ3] : 



dg 1 


11 3 


dfi 




dg2 


1 3 


dfi 


_ 1£ 1 92 ' 


rhi\ 


d ( m 2 


9l) 


dfi V g\ 



L We would like to thank Manuel Drees for explaining this point to us. 
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The solutions of the renormalization group equations are 

9iM 



9i 



rrii 



'1 - 



9i 



A 
A 



mlM + ml(^o) ~ rh 2 + Yl ("^1(^0) - ml) , 

6 

A(fjLo) + 6 (m 2 (/i ) -m 2 ) + — (mi(/i ) - mi) 



5 , 2(/^o)\ 6 / 0i(a*o) 



<72 



#1 



(4.8) 

(4.9) 

(4.10) 
(4.11) 

(4.12) 



where i = 1,2, 61 = 11/87T 2 and fe 2 = 1/87T 2 . Ignoring the running of the gaugino 
masses and gauge couplings, we find that 



Air 2 
3 

3 

~8^ 



m 2 



9l + 



mi 



m 



00 



9l 



m 2\ 2 . / m l \ 2 



Ao [<?2 + 3?] 



(4.13) 



where the subscript denotes the values of parameters at the high scale /i . 

For universal boundary conditions, as in minimal grand unified supergravity, 
the high scale is the GUT scale u x ~ 3 x 10 16 GeV, rhi(fix) — m 2(^x) = ffi 
and gi = a/vt/10 ~ 0.56, g 2 = ~ 0.72. Then we just use RG equations to 

run the coupling constants and masses to the scale of the saddle point /xo = 4>o ~ 
2.6 x 10 14 GeV for M P = 2.4 x 10 18 GeV, = 1 TeV, A = 1. With these values 
we obtain 



K x « -0.017£ 2 , 
AT 2 « -0.0085£, 
K 3 « -0.029 . 



(4.14) 
(4.15) 
(4.16) 



where £ = m/m^ is calculated at the GUT scale. 

Typically the running based on gaugino loops alone results in negative values of 
Ki p6[] . Positive values can be obtained when one includes the Yukawa couplings, 
practically the top Yukawa, but the order of magnitude remains the same. 
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Figure 4: The running of for the LLe inflaton when the saddle point is at <fio = 
2.6 x 10 14 GeV (corresponding to n = 6, = 1 TeV and A = 1). The three curves 
correspond to different values of the ratio of gaugino mass to flat direction mass at the 
GUT scale: f = 2 (dashed), £ = 1 (solid) and f = 0.5 (dash-dot). 

Thus radiative corrections modify a and we need to finetune the potential to a 
few (but not all) orders in perturbation theory. However, although not completely 
disastrous, this can hardly be considered a satisfactory situation, and in the conclu- 
sions we speculate about possible remedies. 

4.3 The inflaton and LHC 

Let us recall that the constraint on the mass of the n = 6 flat direction inflaton in 
Eq. (12T27D reads 

m^o) * 100 GeV ■ A" 1 (j^) ■ (4-17) 

As mentioned earlier, this is the bound on the mass of the flat direction during 
inflation, determined at the scale = <pQ. Since the inflaton mass runs from <fr down 
to the LHC energy scales, it will also get scaled. 



-20- 





1.3 


~~ 1 1 


i i i i 




1.25 








1.2 






*"~o 

& 


1.15 






E 


1.1 


- 


x^ S 

\ S 


eg -©- 

E 


1.05 
1 




\ 

S 

X. \ 

>v \ 
\ N 
\. \ 

S 

\ 

N. \ 
\, \ 




0.95 




\. \ 

\ 

\ V 




0.9 







3 5 7 .9 ,11 13 15 



Figure 5: The same as Fig. 1 but with 4>q = 10 14 GeV (solid) and 4>q = 10 15 GeV (dashed), 
and £ = 1. 

For LLe we see from the solution, Eq. ( |4.10| ), that the flat direction mass only 
gets larger due to the gaugino running, 

m l = m lM + ^2(^0) -^l + ji ~ ™3 • ( 4 - 18 ) 

This has been depicted in Figs. f| and |5|. The x-axis is the scale /1 which varies from 
the saddle point VEV at 4>o down to 1 TeV scale, where LHC is probing new physics. 
We have also assumed unification of gauge couplings at the GUT scale. We find that 
the mass of the inflaton increases only very slightly at the TeV scale. In Fig. |3] we 
see that the increase is \fl$ x 100, x 100, VTJ x 100 GeV, for £ = 2, 1, 0.5, 

respectively. Changing the initial VEV from 10 14 GeV to 10 15 GeV results only in a 
minor modification in the running of m|. This has been depicted in Fig. [|. 

The mass at the TeV scale will increase further if A decreases much below one. 
For A = 10~ 2 , the bound exceeds 10 TeV which escapes LHC limit. On the other 
hand, note that LHC will never probe the flat direction mass directly, but may set 
a limit on the slepton masses. However, we do not claim that LHC can discover 
MSSM inflation, but it can certainly rule out the possibility. If LHC does not find 
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low energy supersymmetry within ~ TeV, then MSSM inflation is effectively ruled 
out. 

The situation would be similar for udd without the top squark. For the Usdd 
direction it is possible that the inflaton mass gets even smaller at the weak scale. 

4.4 A 6 vs. A 3 

One final comment is in order before closing this Section. Unlike m^, there is no 
prospect of measuring the A term, because it is related to the non-renormalizable 
interactions which are suppressed by Mp. However, a knowledge of supersymmetry 
breaking sector and its communication with the observable sector may help to link 
the non-renormalizable A-term under consideration to the renormalizable ones. 

To elucidate this, let us consider the Polonyi model where a general A-term at 
a tree level is given by 

m 3/2 [{a-3)W + (f>{dW/d(f>)), 

with a = 3 — v3 [0] . One then finds a relationship between A-terms corresponding 
to n = 6 and n = 3 superpotential terms, denoted by A 6 and A 3 respectively, at high 
scales: 

One can then use relevant RG equations to relate A 6 which is relevant for inflation, 
to A 3 at the weak scale, which can be constrained and/or measured. In principle 
this can also be done in general, provided that we have sufficient information about 
the supersymmetry breaking sector and its communication with the MSSM sector. 

4.5 Supergravity corrections 

SUGRA corrections often destroy the slow roll predictions of inflationary potentials; 
this is the notorious SUGRA- r] problem [p7fl . In general, the effective potential 
depends on the Kahler potential K as V ~ (^e K<yLp * ^/ m pV{4>) j so that there is a 
generic SUGRA contribution to the flat direction potential of the type 

U(0) = H 2 Mlf (JL\ , (4.20) 

where / is some function (typically a polynomial). Such a contribution usually gives 
rise to a Hubble induced correction to the mass of the flat direction with an unknown 
coefficient, which depends on the nature of the Kahler potential 12 . 

12 If the Kahler potential has a shift symmetry, then at tree level there is no Hubble induced correc- 
tion. However, at one- loop level relatively small Hubble induced corrections can be induced pa, |29|. 
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Let us compare the non-gravitational contribution, Eq. (2^), to that of Hubble 
induced contribution, Eq. (|4.20 ). Writing / ~ (<f)/Mp) p where p > 1 is some power, 
we see that non-gravitational part dominates whenever 

H^Mlfjt-J «mj$, (4.21) 

so that the SUGRA corrections are negligible as long as 0o *C Mp, as is the case 
here (note that H^Mp ~ m^o). The absence of SUGRA corrections is a generic 
property of this model. Note also that although non-trivial Kahler potentials give 
rise to non-canonical kinetic terms of squarks and sleptons, it is a trivial exercise to 
show that at sufficiently low scales, H in { -C m^, and small VEVs, they can be rotated 
to a canonical form without affecting the potential 13 . 



5. End of MSSM inflation 
5.1 Reheating and Thermalization 

After the end of inflation, the flat direction starts rolling towards its global minimum. 
At this stage the dominant term in the scalar potential will be: m^tfi 2 /2. Since the 
frequency of oscillations is ui ~ ~ 10 3 iJi n f, the flat direction oscillates a large 
number of times within the first Hubble time after the end of inflation. Hence the 
effect of expansion is negligible. 

We recall that the curvature of the potential along the angular direction is much 
larger than H? n{ . Therefore, the flat direction has settled at one of the minima along 
the angular direction during inflation from which it cannot be displaced by quantum 
fluctuations. This implies that no torque will be exerted, and hence the flat direction 
motion will be one dimensional, i.e. along the radial direction. 

Flat direction oscillations excite those MSSM degrees of freedom which are cou- 
pled to it. The inflaton, either LLe or udd flat direction, is a linear combination of 
slepton or squark fields. Therefore inflaton has gauge couplings to the gauge/gaugino 
fields and Yukawa couplings to the Higgs/Higgsino fields. As we will see particles with 
a larger couplings are produced more copiously during inflaton oscillations. There- 
fore we focus on the production of gauge fields and gauginos. Keep in mind that 

13 The same reason, i.e. -ff; n f <C also precludes any large Trans-Planckian correction. Any 
such correction would generically go as (7Ji n f/M*) 2 -C 1, where M* is the scale at which one would 
expect Trans-Planckian effects to kick in . 
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the VEV of the MSSM flat direction breaks the gauge symmetry spontaneously, for 
instance udd breaks SU(3)c x U(1)y while LLe breaks SU(2)w x U(1)y, therefore, 
induces a supersymmetry conserving mass ~ g{<f>(t)) to the gauge/gaugino fields in 
a similar way as the Higgs mechanism, where g is a gauge coupling. When the flat 
direction goes to its minimum, (4>{t)) = 0, the gauge symmetry is restored. In this 
respect the origin is a point of enhanced symmetry |J. 

There can be various phases of particle creation in this model, here we briefly 
summarize them below. Let us elucidate the physics, by considering the case when 
LLe flat direction is the inflaton. 



Tachyonic preheating: 

Right after the end of inflation, when we are close to the saddle point, the 
second derivative is negative. One might suspect that this would trigger tachy- 
onic instability in the inflaton fluctuations which will then excite the inflaton 



couplings to matter pTl |33 



However the situation is different in our case. As mentioned, only inflaton 
fluctuations with a physical momentum k < will have a tachyonic instabil- 
ity. Moreover V" < only at field values which are ~ <f)Q. Tachyonic effects 
are therefore expected be negligible since, unlike the case in [plj] , the homo- 
geneous mode has a VEV which is hierarchically larger than (we remind 
that (f) > 10 14 GeV) and oscillates at a frequency to ~ m^. Further note 
fields which are coupled to the inflaton acquire a very large mass ~ h(po from 
the homogeneous piece which suppresses non-perturbative production of their 
quanta at large inflaton VEVs. We conclude that tachyonic effects, although 
genuinely present, do not lead to significant particle production in our case. 

Instant preheating: 

An efficient bout of particle creation occurs when the inflaton crosses the origin, 
which happens twice in every oscillation. The reason is that fields which are 
coupled to the inflaton are massless near the point of enhanced symmetry. 
Mainly electroweak gauge fields and gauginos are then created as they have 
the largest coupling to the flat direction. The production takes place in a short 
interval, At ~ (gm^o) -1 ^ 2 , where 0o ~ 10 14 GeV is the initial amplitude 
of the inflaton oscillation, during which quanta with a physical momentum 

1 /2 

k < (gm^cfio) 1 are produced. The number density of gauge/gaugino degrees 
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of freedom is given by [34 



As the inflaton VEV is rolling back to its maximum value <J)q, the mass of the 
produced quanta g(4>(t)) increases. The gauge and gaugino fields can (pertur- 
batively) decay to the fields which are not coupled to the inflaton, for instance 
to (s)quarks. Note that (s)quarks are not coupled to the flat direction, hence 
they remain massless throughout the oscillations. The total decay rate of the 
gauge/gaugino fields is then given by T = C (g 2 / '487r) g<p, where C ~ 0(10) is 
a numerical factor counting for the multiplicity of final states. 

The decay of the gauge/gauginos become efficient when 

( ^p™^y/' (52) 

Here we have used (<f){t)) ~ (pQfn^t, which is valid when m^t -C 1, and V ~ 
t~~ x , where t represents the time that has elapsed from the moment that the 
inflaton crossed the origin. Note that the decay is very quick compared with 
the frequency of inflaton oscillations, i.e. V ^> m^. It produces relativistic 
(s) quarks with an energy: 

E = I^( t ).(i^) " 2 . (5.3) 

The ratio of energy density in relativistic particles thus produced p re zwith re- 
spect to the total energy density po follows from Eqs. (|5TT|), fl5.3|): 



P —~W- 2 g, (5.4) 
Po 

where we have used C ~ O(10). This implies that a fraction ~ O(l0~ 2 ) of 
the inflaton energy density is transferred into relativistic (s) quarks every time 
that the inflaton passes through the origin. This is so-called instant preheating 



mechanism [35| . It is quite an efficient mechanism in our model as it can 



convert almost all of the energy density in the inflaton into radiation within a 
Hubble time (note that Hr£ ~ lOW 1 ) 15 . 



14 In a favorable condition the flat direction VEV coupled very weakly to the flat direction inflaton 
could also enhance the perturbative decay rate of the inflaton [^6| . 

15 We emphasize that reheating happens quickly due to a flat direction motion which is strictly 
one dimensional in our case. Our case is really exceptional; usually, the flat direction motion is 
restricted to a plane, which precludes preheating all together, for instance see p3|. 
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5.2 Towards thermal equilibrium 

A full thermal equilibrium is reached when a) kinetic and b) chemical equilibrium 
are established. The maximum (hypothetical) temperature attained by the plasma 
would be given by: 

T max ~ V 1 ^ ~ K0 O ) 1/2 > 10 9 GeV . (5.5) 
This temperature may be too high and could lead to thermal overproduction of grav- 



itinos |37, [38|]. However the dominant source of gravitino production in a thermal 
bath is scattering which include an on-shell gluon or gluino leg. In the next subsec- 
tion we describe a natural solution to this problem and show that the final reheat 
temperature is actually well below Eq. (|5.5|), i.e. Tr <C T max . 

One comment is in order before closing this subsection. The gravitinos can 
also be created non-perturbatively during inflaton oscillations, both of the helicity 
±3/2 |39| and helicity ±1/2 states |40| . In models of high scale inflation (i.e. H ini ^> 



7713/2) helicity ±1/2 states can be produced very efficiently (and much more copiously 
than helicity ±3/2 states). At the time of production these states mainly consist of 
the inflatino (inflaton's superpartner). However these fermions also decay in the form 
of inflatino, which is coupled to matter with a strength which is equal to that of the 
inflaton. Therefore, they inevitably decay at a similar rate as that of inflaton, and 



hence pose no threat to primordial nucleosynthesis |4T[ . 

In the present case ~ rn 3 / 2 ^> H m ^. Therefore low energy supersymme- 
try breaking is dominant during inflation, and hence helicity ±1/2 states of the 
gravitino are not related to the inflatino (which is a linear combination of leptons 
or quarks)at any moment of time. As a result helicity ±1/2 and ±3/2 states are 
excited equally, and their abundances are suppressed due to kinematical phase fac- 
tor. Moreover there will be no dangerous gravitino production from perturbative 
decay of the inflaton quanta [^, Q . The reason is that the inflaton is not a gauge 
singlet and has gauge strength couplings to other MSSM fields. This makes the 
inflaton — > inflatino + gravitino decay mode totally irrelevant. 

5.3 Solution to the gravitino problem 

In order to suppress thermal gravitino production it is sufficient to make gluon and 
gluino fields heavy enough such that they are not kinematically accessible to the 
reheated plasma, even if other degrees of freedom reach full equilibrium (for a detailed 
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discussion on thermalization in supersymmetric models and its implications, see || 
H). This suggests a natural solution to the thermal gravitino problem in the case 
of our model. Consider another flat direction with a non-zero VEV, denoted by </?, 
which spontaneously breaks the SU(3)c group. For example, if LLe is the inflaton, 
then udd provides a unique candidate which can simultaneously develop VEV 16 The 
induced mass for gluon/gluino fields will be: 

■m G ~ g(v(t)) < # - (5.6) 

The inequality arises due to the fact that the VEV of (p cannot exceed that of the 
inflaton since its energy density should be subdominant to the inflaton energy 
density. 

If gipo ^> T max the gluon/gluino fields will be too heavy and not kinematically 
accessible to the reheated plasma. Here <po is the VEV of udd at the beginning 
of inflaton oscillations. In a radiation-dominated Universe the Hubble expansion 
redshifts the flat direction VEV as (ip) oc if 3 / 4 , which is a faster rate than the 
change in the temperature T oc H 1 ^ 2 . Once g{ip) — T, gluon/gluino fields come into 
equilibrium with the thermal bath. As pointed out in Refs. H || , if the initial VEV 
of udd is 

<^ >10 10 GeV, (5.7) 

then the temperature at which gluon/gluino become kinematically accessible, i.e. 
g((p) ~ T, is given by [g 17 : 

T R <10 7 GeV. (5.8) 

This is the final reheat temperature at which gluons and gluinos are all in thermal 
equilibrium with the other degrees of freedom. The standard calculation of thermal 
gravitino production via scatterings can then be used for T < Tr. Note however that 
Tr is sufficiently low to avoid thermal overproduction of gravitinos. 

Finally, we also make a comment on the cosmological moduli problem. The 
moduli are generically displaced from their true minimum if their mass is less than 
the expansion rate during inflation. The moduli obtain a mass ~ (D(TeV) from 
supersymmetry breaking. They start oscillating with a large amplitude, possibly as 



16 To develop and maintain such a large VEV, it is not necessary that udd potential has a saddle 
point as well. It can be trapped in a false minimum during inflation, which will then be lifted by 
thermal corrections when the inflaton decays (as discussed in the previous subsection) fl3|| . 

17 Note that the conditions in Eqs. (5.6), (5.7) can be simultaneously satisfied easily. 
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big as Mp, when the Hubble parameter drops below their mass. Since moduli are 
only gravitationally coupled to other fields, their oscillations dominate the Universe 
while they decay very late. The resulting reheat temperature is below MeV, and is 
too low to yield a successful primordial nucleosynthesis. 

However, in our case H in i <C TeV . This implies that quantum fluctuations 
cannot displace the moduli from their true minima during the inflationary epoch 
driven by MSSM flat directions. Moreover, any oscillations of the moduli will be 
exponentially damped during the inflationary epoch. Therefore our model is free 
from the infamous moduli problem. 

6. Conclusion 

The existence of a saddle point in the scalar potential of the udd or LLe MSSM 
flat directions appears, perhaps surprisingly, to provide all the necessary ingredients 
for an observationally realistic model of inflation (Jj. MSSM inflation takes place at 
a low energy scale so that it is naturally free of supergravity and super-Planckian 
effects. The exceptional feature of the model, which sets it apart from conventional 
singlet field inflation models, is the fact that here the inflaton is a gauge invariant 
combination of the squark or slepton fields. As a consequence, the couplings of the 
inflaton to the MSSM matter and gauge fields are known. This makes it possible 
to address the questions of reheating and gravitino production in an unambiguous 
way, as we did in Sect. 5. Since udd and LLe are independently flat, therefore, if 
LLe is the inflaton, the udd direction can also acquire a large VEV simultaneously. 
This gives a large mass to gluons/gluinos which decouples them from the thermal 
bath, and hence suppresses thermal gravitino production. As discussed in Sect. 5, 
non-thermal production of gravitinos is negligible in our model. 

In the MSSM inflation model the mass of the inflaton is not a free parameter but 
is related to the masses of e.g. sleptons, should the LLe direction be the inflaton. 
We have solved the appropriate RG equation equations to relate the inflaton mass 
to the slepton masses at energies accessible to accelerators such as LHC and found 
that LHC can indeed put a constraint on the model: it may not be able to verify it, 
but it certainly can rule it out. 

The model predictions are not modified by supergravity corrections, i.e. the 
observables are insensitive to the nature of Kahler potential. MSSM inflation also il- 
lustrates that it is free from any Trans-Planckian corrections. MSSM inflation retains 
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the successes of thermal production of LSP as a dark matter and the electroweak 
baryogenesis within MSSM. 

The existence of the saddle point requires a fine-tuning of ratio of the soft break- 
ing terms A and m, or the parameter a, as discussed in length in Sects. 3 and 4. We 
dealt both with the case of a local minimum and the case of V'(4>) > 0. We found 
that a fine-tuning of the order of 10~ 9 is sufficient. It is therefore necessary to adjust 
the ratio A/m up to few orders in perturbation theory. However, we find that the 
existence of the saddle point is not sensitive to radiative corrections so that saddle 
point inflation can always be achieved for some value of the ratio A/m. 

However, it is conceivable that the mechanism of supersymmetry breaking, which 
lies outside the effective theory of MSSM combined with gravity discussed in this 
paper, could remove the fine-tuning in some natural, dynamical way. For instance, 
A/m could turn out to be a renormalization group fixed point so that once the ratio 
is fixed, it would remain fixed at all orders (for example, see |@). This requires a 
detailed investigation, but it is warranted by the simplicity and the apparent success 
of MSSM flat direction inflation, which is unique in being both a successful model of 
inflation and at the same time having a concrete and real connection to physics that 
can be observed in earth bound laboratories. 
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